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This paper is ccncerned wioh one rwc-dinens icnal supersonic 

flow of 2 ohick ourbulenu 0001 : 030 "/ layer over a oraun of relaoively 

sr^ll wave- like protuberances . The flow ccndioicns and one 

geometry are such ohat ohere exosos 2 sorong inoeracooon beov/een 

one viscous and :.n viscid flow. The problem cannoo be solved wioh- 

out inclusion of intaracoion effects due 00 one occurrence of 

the separation singularity in classical bo’cndary layer r.eohcds . 

Here the interacting boundary layer equaoions are solved 

n'noeriCcLlly ‘using a tire- like reiaxaoion meohcd vioh ourbulence 

effects reoresenoed by the incl’usion of the eddy viscosity 

nodal of Cebeci and Smioh. Resulcs are presenoed for flew over 

a train of ‘up to six waves for Mach n’umbers cf 2.5 and 3.5, 

Reynclds n'jmoers of 13 and 32 x 10®/neter, and wall tamperacure 

raoions (T /T ) of 0.4 and 0.3. Lirioed ccmoariscns wich 
w' 0 

indapendeno e:cpeririental and analytical results are al 
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ABSTRACT 

This ?ac«r is ccncarr.ad wzzr. the tvo- 
diaenaional iucersar.ic flsw ar a thiajt tur- 
iulant: aounGar-/' layer aver a traxa at reia- 
tivaiy small wave- lake arrfrherances . The 
flaw ror-ciGians and the teomacry are such 
tmar thare exists a strcng ar.teracticr. be- 
tween tea viacous aad i.tvtscid flow, rhe 
orebiam rar*r.oi: be solved without inclusion of 
izitaraction ef facts lue to toe occurrence of 
tne seoaration sincularitv in classical 
ooundar'/ layer tethods . Here the interacting 
'couadar'/ Laver ecuaticns are sol'/ed numeri- 
cally isinc" a time-Lixa relaxation tathed 
with" turbulenca affects represented oy the 
inclusion of the addy •/isccsitv tocal of 
laceci and Smith. Results are presented for 
flow over a trai.t of up to six vavas for 
Mach fomDars of 2.5 and 3.5, Reynolds n’umhers 
of 13 and 32 x 10° -tetar, and wail tempera- 
ture rations '.2 /T ) of 0.4 and 0.3. limited 
V o 

ocmpariscns with independent experir.ental and 
analytical results are a«so piven. 

;ic.ME:iciArr?.£ 

a .Amplitude. 

.A iddy viscosity iampinp f'lnction. 

Ga skin friction ooeffioiant '.ya* -*“/2. 

G Constant oressura soeciiio neat. 

0 

r Momalicad longitudinal valocity, 

7 '.:orma_icad total enthalpy, g » S/H^ . 

.1 Seat transfer ooef f ioier.t . 

2 Mondoiaensior.al total enthalpy, /u^' 

K, Constants in eddy •Ttscosity models, 

i '/iscesity parameter, . 

I Mixing lengt.h, 

1* Raferer.oa le.tgth. 

M .Mach n’lmcer. 

o Mondi.mens icnal static oressura, 

0 - o’/o* u'2 . 

» » 

?r Prar.d^l n’JurJser. 

?r,„ T'-ixbuI.ar. w r*’JULbiar. 

L 

c_ Mor.di.tensional turbulent heat flux 


Reyncids number based on reference 
viscosity, Re„ • ?-a, u^/u 1 (u*2 /c*) , 

Ravnolds n'lmber oasad on free stream 
• _ * * , » • 
viscositv, ?e a J U, - , . 


Mondimene icnal static temcerature , 

2 > 2'G'/u*2 . 

3^ * 

Mcndimensional x. and velocity 
oompenents , i»u’/u*, v^v* ReX'’'u^. 

2 rans formed v velocity in tne ooundar*/- 
lay er . 

Mondimensional coordinates (surface or 
Cartesian) , £, = P.e-' * 

Wavelength . 

1 

u;/ 2 ^ ■ 

Pressure gradient oarametar, (22/u ) 
(du^/ds) . 

Ratio of specific heats , 

Transverse .ntemittancy f'cnctior. . 
.Mondi.me.t 3 ional displacement thic.tr.ess . 
-ncompressibie dispidca-ment thicxnass . 
Oispiaca.te.nt cody haigr.t. 
iddy viscosity. 

eddy viscosity pararu*.tar, - ' . 

. o- 

Eddy viscosity parameter. s^L- y 5 -^ 
Transformed .tomal variaale. 

Static temcerature ratio, -s => T/2_^. 

Surface inclination of the iispiaca- 
mant body . 

longitudinal inta.rmittanoy f’lnction. 
Mondimensionai viscosity, 

.“•j */'j * , u’ 2 /,-*5 . 

Transformed longitudi.nal variaola. 

r'lnctionai grouping l.n i.tner recior. 
eddy viscosity model. 

i’ondi.mens ional turbulent shear stress. 

• « 

Mondi.tens:.onal iensitv, j » 3 o _ . 
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e 

Conditions evaluated 
body or at me outer 
ary layer. 

an the displacement 
edge of th# bound- 

i. 3 . 

Flat plate value. 


1 

Index for tne longitudinal finite dif- 
ference nean. 

w 

Conditions evaluated 

at the wall. 

» 

Conditions evaluated 
f reastraam. 

in the upstream 

3a=e 

rscriots 



* l^acocaa iiisaiiaianal qua^iricies . 
3THCCUCTICN 


- • Sour.dar^/ Laver £juatior.s .n Phyaisal. 

Cobrdi.iaras 

The suiraJaility of the interacting 
soundarv layer ecuacions far describing the 
relatively strcng streair-ise variaciar.s in 
*^e boundary *ayer characteristics due to 
sudden chances in the ccdy gecmetrt’ has been, 
least for me laminar case, veriiiad 
'^seller (1, . This looroach is ised .n 

the present study in vnicn ?randtl ' s classi- 
cal boundary layer equations ere adootad 
VI th the only modification that me pressure 
’^^riaticn was not prescribed but calcu-atad 
simuitaneously from a viaccus-inviscid 
interaction ncdel. 


This pacer is concerned with the two- 
iinensional sucarsonic flow of miok t’lrou- 
lant ocuhdary .ayers over a tram of rela— 
cively smai— wave~liJ<a oroturarances . In- 
terest in mis subject arises from the need 
to prediot me extant to wnioh an initiailv 

plate boundary layer has been aisturned 
by a regular corrugation in me wall surface. 
The flew conditions and the geometry con- 
sidered nare are auen that mare a.-cista a 
strong interaction cetveen me viscous and 
i.T’/iscid flow. The proclam cannot be solved 
wimout including i.nteraotion affects because 
classical boui.dsry layer .tetnods would tar- 
^ i separation pcint singuiaritv. 

To handle t,he present subject by boun- 

netheds, a tacnnicua for treatment 
3f Che i.ntaracting boundary Layer equations 
well as tocels for turoulanca ana for tne 
•/hscous-inviscid interaction process m.’ust 
oe availaoie. A .Tumarical method for addres- 
sing closed oucbla secaratior. regions was 
ieveloped by Verla and ’/atsa (1) . Zt was 
appiuad to a number of Lami.nar~secarated 
flow prooiams including flow over a train of 
sine-wave profioarancas (2). This .tetned 
uses the intaracting bouncary Layer equations 
with a time-iixe rei.axation concept wnior. 
accounts for cha boundary-value nature of the 
proola.m. This approach is adcotac in the 
present study wim. ths i.tci’usion of me 
eddy viscosity modei of Caoeci and Smith into 
the solution scha-ms . Ths present form of ma 
.numerical algorithun includes several modi- 
fications to mac of tne eariiar woric (2, 3) 
in order to acconr.ocace the turculant nature 
of the rlow, the thioic boundary Layer, and 
the rather dramatic geom,atry variations of 
t.te wavy wall . 

2t was fo'md that the metnod was oaoaole 
of handling t.he i.nteracting turbulent flows 
o. present interest. Solutions were obtained 
for flow of thir!< turbulent bcundar'/ layers 
over a trai.t or waves. The results are ore— 
sentad i.n tarms of surface pressure, ski.n 
friction and heat transfer distributions. 

The predictad trends are compared with avail- 
able analytical results based on small lis- 
f’rbance t.heory and with e.'cterimental data. 


^ The boundary Layer approximation i.n tve- 
iima ns ion a., viacous flow prooLams impLias 
— ’'•® pressure variation is assumec to 
occur only along ona cocrcinace, ta.'cen in me 
general direction of the wall shear Lavar. 

The degree of tnis approximation daoer.ds on 
choice of me coordinate system'. vhi..a 
for very mm boundary Lavers over a corru- 
gated wa.,.1 , or thith bcundary _ayers over a 
..elatively clat wa*i, surface coerdmatas 
were suitaole, (sea .=.ef. 3) for chioX bound- 
ary Lavers flowing over a'smali amtiitude 
wa_» , Cartesian coordinates were founc 
to oe more aptrcpriaca. Acoordinciv, the 
governing equations will first oe wrifTen 

oc th tne usual surrace coerdmatas 
>3*, n*) and me 2artasian coordi.natas x* , 

V-) uamg the notation x^) m ianote 

aithar of these. ^^oncilmen3lonai variaoLes of 
order one are now cafined according to ma 
scheme. 


X, 
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ana u , v , p , j and T represent the mean 
■'sl-cities , pressure, der.sitv and tamtera- 
tura respectively. 


The turbula.nt boundary Layer aquations 
im these variacies are: 


r.uitv icuation 




( j v) 


■Mcrr.e.ntum louation 



du 

)u 

® 5 

dx, 7^7 - 

3X, “ c 


C) 


j.r.erT/ Sruacur. 

i'T l'T> 

3 ( u * V 

5X, >X- 


- 3 U — 

a a tLx, 


whera and are the r:or.d‘..~nansicnal rur- 

buienc arreas and ruraulanr beat ilax respec- 
tive iy. 

the qas is assunad tc be air wim ron- 
snanr apeciiic baacs and rcnsranr Prandti 
-LSfflar, ?r » 3 . t2 virb tna perracr jas law, 

State trcaricn 

1 

P * — ^ at (51 


Poundar*^ 

'Oendi 

tl 

ens 

■1 ( .*c. 
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v(x. 
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and 

u(.x.^ 

, ;<,) 

9 

■4,(X, ) 

at X., 

7(x, 


9 

r^(.x,^) 


equation. tn an effort to better align tne 
predictions for separated flows wit.b exoeri- 
r-antal fata, previous investigators (see 
?.efs. £ and 3_ for exa-tples) nave e-tpirically 
Todiriad tne equilibrium eddy viscosity 
nodal to account for t.te nistory effect, 
thus ' rroce.". ' , ‘ reiaxatio.n ' , and ocher ncdels 
were devised and success fully applied i.t 
several of separated flow predictions . tne 
of tne present autnors (o)'also used tna 
'frozen' and 'relaxation''” tocela i.t t.be 
interacting boundary Layer ecuatior.s for 
separated rlows with no signifioant improve— 
in the predicted results over tncse 
obtained with the bas.’.o eddy-visccsitv todei. 
-- appears tnat to achieve tore satisf actor*.’ 
results, a turbulence .tcce_ emplovi.tg the 
P--3*ilent transport partia- differential 
■equations will nave to oe developed for use 
in strongly int.eracting and sacarated 
boundarv layer oalo’uiations . this will, of 
o-ourse , f’lrther ta.x the computing times 
required for t.tasa talc'ulaoions . 

amove tentioned limitations in 
tind, the basic form of the Ceoeci-omith 
eddy ’/isccsioy .todei ’was adaoted for tne 
present stud*/ "wnere i.ttaraotion effects tend 
to reduce longit’udinal gradients and oniv 
^*^i— aecaration regions are encountered, 
thus 'we taxe 


■where xi (x, ) describes the bod*/ surface con- 

tour (x-, =• 3 i.n surface ooordi.tates , xi » 

■*w -V 

tartasian ocordinatas) . 


2 . t’ urbule.tce .'■lede- 

■"0 octain oios^ure of tna 5*/stam of 
equations (2 -t) , todais for the t’crnuia.nt 
stress anc furbuiant neat fluic ‘.erms are 
needed. the add*/ -/isccsity ooncsoc used in 
con^'unction wiih ?randti ' s' tixing length 
n*/pothesis far the wall layer region is the 
tost -widsi^/ used algebraic node! f-oc fornu- 
-®ut stress, .i well 'tnown representation is 
the fwo layer edd*/ viscosity toda- of tebeci 
and icith wr.io.b has been ver*/ sucoassful i.n 
.todaimg t'urnulence effects far flat data 
ooundar*/ layers end otnar attached boundar',* 
layers ■with moderate pressura gtedients. 
less ravoraoie results are obtainad wr.en 
-‘Si.ng mis todau for strongi'/ intaractin— 
and separated flow regio.*us wnera it appears 
to fail oor.captualiy . 

2n general turbulent gua.ntities liJ<a me 
Reynolds stress are governed 'ey transport 
equations, tn'is requiring that me furbuler.ee 
.nistory be accounted for. The add*/ ■/isoosit*/ 
concept reiatas the Pmynoics stress to onl^/ 
t.ba local r.aan flow gradient. This oorre-* 
spends to t.te physical idea that production 
of furb'uience at a poi.nt i'ua to interaction 
■wim the mean flow is oar.cailed by the dis- 
5i?d-ior! due to its self-i.nteraction (mis 
13 referred to as the 'local equilibri'um" 
concept) . r.n ot.ber ■words, t.be eddy viscosif/ 
model 13 t.he solution to a tr*u-"'.cated transoo'rt 


and relate c„, to t_ oy furbuient Prar.dtl 
numaer as ‘ " 

■ *T ' "r jx, ’ ' ^-•■r ix, ' 

fursuie.nt Prandtl n'umner is nera ta,<a 
constant, Pr,, » 3.30. The fwo .aver (out 

and inner region) Ceosci-^amith model is : 
given as : 


r.n.ner P.ecion 
(e/u) = 


’ -'2 I *1 

i I j'u ! 


■wnere i 3 .K, x,,;'! - e.-cp (. ); 
■wim :<, *3.40 and 


• • • ♦ 1 !1 W 

A > 25 .U /u ) (u , , ) 

, 3X., !*w 


■Where me absolute value of j'u nas 

been introduced in ec^uation 3c) as "a mod: 
4ioation of the Ceoeci-imim mode, for 
reverse flows . 

Outer ?5cior. 

♦ w 

( € / 'J ) . = — K . V • , 


REPRODUCIBilunT OK I'liK 
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wnare v is transverse intemittancy 

function 

; - (I - arfC5(x./x^ - 0.73) ]}/2 (9b) 

The veriaaia is t-ta value of x, at whicb 

'u/u ■ 0.995 , and s. . ^ is the incororessibla 
a xinc 

dissiacament thioicnass . 

3 . Soundar*/ laver iruations in Transforr.ed 
Varianlas 


where 4 is t.he viscosity oaratetar defined by 
i - 3u/3^4^ (15) 

with 'j given from Sut.tarla.tds viscosity law 
and the turhuiant parametars , 2 and i are 
ia fined as 


c » 1 *• (t/u) r 


£ » 1 - (c/a) r 


(15a) 
( 15b ) 


The boundary layer aquations given m 
Section 1 are here recast using tna lev"/- 
lees transfomation. 

The new independant variaoiae are ie- 
finad by 


wr.are .' is tha atreamwisa intamittancy 
f'unction : for full/ lamnar flow * » 0 and 

far fully tuaouiant flaw ■ 1, wnila far the 
transitional region its va_ua varies amcotr.ly 
from zero to one. The parametars i and i 
are obtained from tha local in-'visaid flow as 


0 * ■ '-10 


0 ix. 


( 10a, h) 


The normalioad fasendant variaoies are now 
defined as: 

'/eiacity ratio, 

? » u/u, (11a) 


a . u^-T^ 


-- 


State to'iation 




(17a) 

(17b) 

{ 13a) 


near, static a.ntnaiov ratio 




3 , -3 a -2, 

— r (3 - -r- .- / 


(13b) 


(llh) 


!r, mean total ant.talpy 
? - H/2. 


(lie) 


Wit-t these definitio.ns , equations (2-4) 
become : 

Continuitv tcuatior. 


3V 5r - . 

— -2? r-r - 0 

3 n 3 % 


(12) 


Momanfum 3cuation 


Id. - 7 

) n 


3 (3-?*^) - I- U2 I4-) 

in 


or 




il 

3 n 


(It r) 3(g-?“)- ^) 

a an j n 


(12a) 
3 3? 

a.”" ■ 

3 n 

(13b) 


atatio Ta.trgerafure E.nerg*^ louation 

:r r 4 . 7 i± . , 1 ^ ^ (i 11 ) 

in 3n in ?r)n 

or (1^4) 

Total Temoeratura inarf/ 3cuaticn 


2C r ^ ^ 7 


3n 2-; 


(4(r - ;/2.r): 




- (ii. ii) 

3n '?r 3n 


(14 b) 


Boundar’/ Conditions 
.=• ( i , .0 ) *0 

7(:,n) =• 0 

9 ( i , T ) 


(13a) 


* -v^ 

or g(-;,.i) - 3,^/3, 

where » 0 for surface cocrdinatas and 

ar 

* '..(5) Cartesian coordi.natas . Also 

V V 

we have t.hat 

.=•(-:, .1 ) * 1 

3(i,n) =1 as--- V 13b) 

or ? (< . T ) - 1 

The turbulence ra-atior.s given in Sec- 
tion 2 can be expressed i.n transiorrr.ed 
van sales as: 


2 2 2 2 2 

, , , — -a -'-1 ■’ ,1? 

(c/u) . = >?e^ ■ ~ 1 ~ (20a) 

— r , , 3 ' > ■- 


■J ' - 1 4 3 
3 


Where 


1 - ixp(~r^) 


'20b) 



X, 3^ U 
- i i . 

^3 ^ V, 

:6i 3 - 3 ^ 
Ou~ar .^acion 

(s/u)_ > -2 — 1 


3 « ?.e 


■•■. ‘v ,!F ,--'2 

7 ;;^ tt: > 


(20c) 


r 5 . 


?.a_ 


<inc 


i. 


<iac 


?a _ 3 j 
r a 4 


i 3' 


3C-?) dn 


(20d) 


( 20e> 


Nora mac, aa far as me fsm af ma 7 cverr.- 
a<7-ac3.sr.3 .3 acncarnad, ma aniy dif- 
faranca aec-/aen m.a aae ar surfica ccordi- 
.racas and >.arrssian acorddcacfts is m me 

3«'’indary aor.diraan aquarian 13a) -m-.g 

can aa alinanarad laint? Jranar.-a tranacosi- 
cian maoram ay wrizing mac 


7 . 7 - 2 ^ T ’ ? 

w 

Vim maaa cransda rrncian acuacians 
ind (13) ’/laid; 


ar 7(1,-) - 1 


(22b) 


cncaracci.cg acundar*^ layer 'alc’la- 
caons rsquira an initial /eiicity and tar- 
peracura arafile at aotce scacian ahead ai the 

rsqian (see riq-are 1). 
...13 ?rwci_e was abcained ners fran a non- 
-.ncaracci.ng rwo dirte.nsianal laainar-transi- 
Ciorax-curbuienc 'aoundar-/ .aver aalcula-'a-' 

a^es"^-^’";-"^-”' ^^--=•‘'•"•■■9 tecmiqua is^g''! 
?.escr-red arassure diatribucian . 


(21a) 
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(21b) 


(210 


14 ) 



■U/'J.. =1 


12 ! 


j/u. 


1 M/ 3 - 2.53 

1^602-10.32 X ,0?'m 

=0.31 


Initial VeicoTy 
-^rof iia 


Canci.tuitv ^cuacian 

il ~ 2^ ^ ' 

SiT " )7 
Mcntenc-m truacian 


(22a) 


or 


T 7 i 


7 

3 ' 


7 - 3(3-?^ 

^ n 


1 - (is -^-l (22b) 

3 n in 


7 il 
3 n 


xnerc-/ rruacian 

■ r hi .. y hi 
3 i 3 n 


( 1 - ▼) 3 ( 7 -?-)- ill 

' " 5 n 3 n 

(22c) 


zi2 (i2) ■- i_ iii 
3 n 


3 n 


3n (22d) 


ar 


)l ? h£ - 7 i£ 
3 i 3 n 


U(r-s/?r)- 


, ih 
3 n 


3 1 ■ ~ 3 1 


Boundar-/ Canditiar.s 


(22e) 



4 ^ -^1 .. 

---s* I .-lew jigcr.err*/ 


ZTj.Z'B T 3.C* i- (TT- '*^cce ' 

The i-ncaraccian af -me aou.ndar'-’ lav-- 
vim ma isentracia aurerscr.ia mviscid' flaw 
la ccca.ac in t.-.a tressura grad.ant ta-a-;--ar 
3 by aauplmc it ta me ..tci.nacian C 
c.ne tacal diaplnca.r.enc aod',- Sm. Th® 4--- 
prassura 13 aecainad fron tna ' Pranatl-yever 
-2-.iC-.cn ipcrcxiciac2G .i2ca oc s^^cr*^ 
in -srr.s cf 3-, as ” " 
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vhera - can"' (do,/ ix) 


and 


? ( j , 0 ) 

3 

'^ ( i , C ) - D 

tf 

3 (C, 0 ) 

m 

3 ,(i) 

0 ■» 

7 (L 0 ) 

9 

7 .^ ( ; ) 

Cla) 

Once a. 13 ■ 

? ( i , ») 

9 

t 

ind iuiar's 
ecuati an , 1 ' 

3 ( i , » ) 

9 

1 

VI 3 ecus fla^ 




:a isentra: 

L r ® • « 


3i 24) 

( 25 a; 
( 25 b ) 

'250 

ralatiar .3 


since they ire diractlv aar.n-ctad' 


•3 



(27) 


diiolace.'Mnr thickness given in equation 
(2Sc) . 

N’w'MSaiCAi, ;-<STHCD OF SOLUTION 

The n'laarioal .tethod isad is an inpiioit 
finite difference schema written for the 
similarity fora of the governing equations 
that marches from some initial station alona 
the surface to the terainal point of inter- 
est. To account for t.te boundary value 
nature of the prco.ea, Verie a-td" Vatsa (1) 
have added t.he tire dependent concept, “ 

3 i m lar to tne one ised for the solution of 
a_liptic partial differential aquations. 

This results in todificaticn of only t.-.a 
momentus aquation .lib) oy repiaci.ts the 
pressure pararetar i with"! defined as 


1 -1 


I . , 


(2S) 


This methcd has cae.n successfully aopliaci to 
laminar separated flow prcoiems with 'arious 

concigurations includinc one with mul- 
tiple .nteracting regions (2,' 3) . The 
e.xtension of t.his apprcac.t to turoulent 
boundary layers i.tvoives, aside ftom inclu- 
sion of the eddy -/xsccsity medal into t.te 
solution senema, a numaer of modif ications 
(see auso ?jzz. 5) . Specif icaily , the 
-Cllcwing steps vera taxan. 

1. The r.imiarioal staoility and conver- 
gence rata nas beer, ennancad bv i.ttroduci.tg 
a new diffarencing .n tne ocntinuity aqua-’ 
tion. It nas only recently been recogr.icad 
tl' --e ’-cncitudinal darivatives in tne 

ocntinuity aquation provtda a pacn for 
intaracting flows to propagate' i.tforaation 
7o icoommodaca thi.s -.n-of- -•» l_y 
req’nras the use of some sort of a forward 
i-ffarance procedure. *.n tne present worn 
we acopt i.n the conti.tuitv acuaticn the 
fcilowi.tc forward diff arer.ci.-.g 


c) 


• 3F, 
.v-r) 
3 1 


whare tne superscript ^c) danotas values at 
® pc'a'rious time stap anc subsenot i refar; 
to the _th station alo.ng t.ta langt.t of the 
surface . 

2. The reliaoility of tha peasant 
alcorit.tn was a.thancad by adcoti.ng t.te ' uo- 
wind du.f fa rare mg' ocncapt for tne lor.ci- 

-snvecticn artacts. In tna reversed 
flow region upwind diffarencing was usaa 
t.te Icngitudinal direction for the convective 
tarms in orcar to satisfy the staoility 
requirs.me.tts. This eiiminatas tne sb-oallad 
artioicial convection' ocncapt used eariiar 
(2) for the -aminar case. This .mcdif icacicn 
3igr.ifica.tt oecausa tha velccituas i.t t.ha 
reversed flow regions are larger in t.te 
turbule.to case tha.t i.t t.ta lami.tar. Thus 
-..e convective term i.t t.ta .mcmentum aquation 
IS diffarenced as 


> 1 


> (T, - F...)/ui - 


V,|) 1?':! - 


IS replaced bv 


(o) 
■ 1 


forward flow, and by 
3y rapiaci.tg t.he 


for reversed flow. 

with ?/'** t.ha occurrence of a separation 
point 3i.tgu..arity is avoided ( 1 , ' 5 ) . Mote 
t.hat tha first tara on t.te rig.tt tand sida of 
aquation (2-) vanishes ror reversed flow, arc 
^a second tara vanishes for tne forward 
flow. The same procedure was foilewed with 
t.ta tara Fj-J/n m. the energy aquation. 

* — — tr Q rto re , tne uowmd differsncin— 
was found also helpril i.t the -i direction, 

convective dcmi.tatad outar ragiots of 
me t.hic.< turoulent bcundar- laver. It was 
oo our ittuntion (3) t.hat upwmd 

? 3f the j?'»n tara i.t tr. a temen- 

the cenver^anca mzarta oz the .t'^erical 
(see also .^ef. }_) . I.t the ooundar.’ 
layer near me vail t.ta diffusion tara .' ' 
dominatas the oor.vectu ve- like term ?. and^ 
a central diffarenca sc.ta.me for .s 

.icwflver*, the tuta.? raac*“*s 
oz zz,& boundazy Ijv^r vie di,di'j 3 un "s*^ 
decreases 3 LZTi. 1 . f \ tan" J, v inci .'I’vnei"' ■“a.l *"• — 
stoDi-iCitis tcc’ar. Trcti i snudy ad -he Ticcai 
aquation ^ » 3 _t is found tnat wit.h 

central diffarencing t.te oritaria i m | <1 
must =e adhered to, to avoid tnase osci_lJ- 
tions . rianoa the term. was rar.oral iif- 
1 in ■ <1 jiid uowi.tc did - 


rarenced when 
farenced when 


3. Tha oonvergenca rata of tna ti.ma 
-Suaxatioh so.i.ution metnoa for tne tnio;t 


ooundar/ .avers .tas oee.t found mu 


f 2 r thi.t ooundart’ layer: 
iiffar largely in tnat fo; 


slower 


Ihe two oases 
:a mio.k oou.tcar-.' 


J ■» m . 

.aver me di.jturoanoa oo the total tiso.aca — 
cody trom tne olat olata val ua * s / 
sma.k_. .t was arguad t.tat tna r.umerioal 
trancatio.t error oan oe of ma sa-o.e order as 
relaoive o.har.ge oer one itaratior. , *.-us 
•reading to '*ery smal. oonverganoa rata. w.» 
mtroducec tnarefore a tew var- - 

place or tr.a total iisplacament oedv 'jm. 

Th . Dm 13 tafinad as T_ a [ a .. ( i , t,' -i - C " ' , 

whera us tne iisplaca.t ent mic.kness 

•-■'•ioiai station a-td .s a constant 
■ 3 rdar of t.te ta.xi.mum profurera.tca 
laloulatior.s performed -iizr. t.tis 
icaticn shew imorcve.mer.t ..t the ocr.vp*- — 
ga.tce rata. 

p-.a accuracy of -he oalcu.atad scluticns 
^6p^nds an th(3 ascirss ad pr^cisi^an 
-i.tita-dirreratca apprcxi.maticr. a.td tne stap 
- ' • — " turo'Uiar.t councary layers laroe 

0. hanges occur i.t t.ta .'“loci tv orofila i.t tha 

1. t.ter layer ver*/ tear t.ta surface. i. su*'®' - 
ciant r.'c.moer of r.es.t pci.tts are needed tear 
t.te wall in ordar to 'get a good oasolction" 

-h _..e .redictior.s oo vaii shear and surface 

**^“^*“** tr.^i auasr ad -r.s 

^C’mdary .ayer wnera the lamz-lees variacla 
' acuires large values, t.ta 'o.ta.toes ara,~on 

■ar na.td, very 3m.-il. This is asoeciallv 


the ot.t: 





in ihe 3asa af a miCK vurbuiana 
aoundar*/ layer -tiarumad by a raiaciveiy 
amail prariiberanca . Thus, far reasons af 
«iii-iancy and accuracy a variable .r.asr. siaa 
in tne i direction la used in aoiving nest 
furbulent coundar/ layers. A nesti proving 
in aiae fro.n tbe vail as a gecnatric aro- 
gressian is ised in tbe present algorit.dn. 
31attnar (^) has s.down tnat in terns of a 
trar.sforoiaa nomai variaaia M(.t) reciacing 
t-da stretebed lew-lees variaoia -i , the 
tr'uacatian error is proportional to as 
- 3 , ar tbe netbod af ' aalculatian is second 
ardar accurata. At tba 30 b grid point tna 
physical coordinata is aotainad frea 


(X 




•L)/(S 


1,'lN. 


•L) 


(23) 


The governing equations ware lineariaad 
and the partis, derivatives ware replacad by 

The addv viscosity taros 


L.tita 


-li «Va-, 


rtarences . 


•:/ a , 


appearing as a nonlinear tarn 


governing equations . is approxinatad by its 
previous station value. lantrai diifarencas 
were isad to represent parti ils vltn respect 
to n (axcapt as noted acove where upwind 
iid^arsncing for ?- was required in t.te autar 
region af the ocundary layer) as veil as far 

the tarns d'i^^ dxr and far d5w --x, , af tba 

pressure gradient cal aulstians . 2 cwmd dif- 
farencing was ised in tna ror.’/ecti’/e tarns 
in tba nonar.tun and energy ecuaticns and 
^ard diiiarancinc 
rontmuity eauaticn. 


where X» -(j-L)iM, (;_ 

e « ^ J *T13X 

» 1, ^d where ill is the constant stap in tba 
“ransf om.ed plana. The second ardar accuraev 
13 acb.aved ay varying and .lolding iM. 

fixed 1_^) . It was found .tare that id in— "" 


staad one raclacad i;i. 


IM, where iri 13 of 


course varying vitn ir-ax' ’•'Ri-a .-.oiding ;< 
il-«ad tna error iimnisnes witn iri as 

^ , i.a. auen faster tnan lil . Figure 2 

snows tba surface rsaating paranetar's dape.n- 
danca an Hi and thus provides an accurata 
error esttnaticn procadure. 3asad an tnis 
step Sica study it was found that with values 
3= 200, 55 , and .X » l.:S4, a 

■~ur.cation error was incurred in the aalcuis- 
tian af wall neat transfar. This represents 
an accactaole tar.promise aeewean the accurar/ 


and the effiaientr.* af talcuiaticns . 



IN^ X iC* 


eg ^ 
.317 y 



Fig. 2 Accuracy Study af ^'/all Heating level 


The ta.cu.atian rcntr.er.cas with certain 
initial ccnditians and tbe.n through tne ti.te 
iependant approacn (1) tne steady state scii- 
tion for a given set~af boundary tonditians 
13 sought. In t.ba present calculations tne 
loitiai conditions vere set ay taxi.tg t.te 
laro tir.e iisplacerant aedy ta correspond to 
i flat piata bouncar’/ layer and tha surraca 
protuneranca ta ae af tero neignt. Sunsaquar.t 
tir.e sweeps are conducted wicn tne wave 

air. plltuce increasing gracua.ly by a sna.l 
ar.ount. Adtar the cesired gaoretry is raacnad 
(after the first 10-15 sweeps) tne" tire- Lixa 
relaxation process is tcnti.tued 'intil t.te 
:lcw properties are reiaxad ta tnair fi-ta^. 
value. This process .s snewn in Figure 3a 
where the sxm iunctian coefficient at one 
location (s 3.53) _s snewr. as a c’uncticn cf 

ti. te Iteration curccer. This lacatian is tear 
the :’inctian af tha flat plats with a singia 
31-te-wave protuneranca vnere separation 
occurs. The resuiti-tg sxi.t crictian and 
surface heating distributions are snewr. i.t 
Xiq'ures 3b and 3c respectively. For this 
case wit.t a thi.t coundary layer, t.ta calcula- 
tion was cerforr.ed i.t surface ccom.tates . 
j"iTure 3a shows that cnce the cull atofuner- 
ance haigr.t is attained .11 sweeps) it tj-xes 
about 50 -tore sweeps for the sxi.t friction ta 

Its 'steady stata' ■.•alue. This ca_- 
*'^»at.cn , wi .fi ^ ^ to ma ^ m a ao m 1 3 and . 
langitudi-tal grid aci.tts was tarfarr.ad in 
5 tinutes af cenputar ti.te an t.te 3T5-L53. 

FJISULrS .Al'ID DIuCrSSIC.N 

•A rsijor i.tterest of tna present ..tvesti- 
gatiar. is i.t t.te turerica. precictians far 
t.tic.k turbulent aoundar*/ layers aver a waw 
wall, as these i.t t.he e.xperments cf Pefara.tce 
( 11) . The gecr.etrv and t.he flew conditions 
were therefore c.tcse.t to coi.tcida witn t.tase 
given i.t P.eferenca (11) . Tha utolituce end 
wave langt.t irs e* ~J.23 cr. , w' * 3.55 cr. 
respectively. .A reference length 1' * 15.25 
cr. was c.tcsen. The base flaw conditions ire 
defined by =• 2.53, ?,e = 10,32 x 10^/t, 

T * 17i’X a.td T -• T * 0*31. ’“.'•ca f "'r-.t , 

we refer ta t.tese canditians as standard 
flew conditions. 

Ta abtau.t the present results it 13 
first .teCQssary ta re.terata ..titiai prafiles 
at so.te ocir.t ir;aad af the first aroficerance- 





■^1^* 2.5 

^6(33 =1.55 X 10 ’’/m 
=0.31 
S = 3.53 


3i4n coordi.nares version oi ziia ocar.dar*.^ Layer 
equations .tore rsasonaoie. The results* of 
t.te calculations sncwn here vers oeriorr.ed 

a Longitudi.nal stepsize ix ■ 0.02/ and 
a 35 point jrid across t.he ooundarv Layer. 

Zxamoies from tne calculated results 
are presented for flow over a train of uo to 
SIX waves, for y.ac.h .■'.•cnoers y. "2.5 and’ 

3.5, for Reynolds n’ntoers .Re "» 10.32 x 10®. m 
and 32.46 x 10® 'm and for wall to total 
temperature ratios T /t, - 0 . .;o a.nd 0.31. 
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'iq. 3 Surface Properties for One Protuaeran; 

(a) .'lumarical Cor.’.’ergenca ?.ata 

(b) Surface Heat Transfer 

(c) Surface Skin Priction 


j '-1*4C*U^3 . FOT -«Q 

conditions tnis station was taken at 
'2.90, where the initial profiles we: 
tai.-.ed from a r.oni.nteractl.-.g calcuia: 
-screspond to the boundar*/ laver as • 
ops along tne wall of t.he* ITWT langla 
Tun.-.el (^1 . The i.-.teracting algori; 
subsequently a.tployec he tween this i.t 
station and a downs trea.-n station oast 
.ast orofuberar.ee . The oroble.t wks 
-omuiated and solved in the rustomar 
-see coordi.tates . It turned out that 
geometry extremes taka tne use of the 


:n to 
de ve i - 
Wind 
r. was 


?ig. 4 3isoiaca.tant 3ody and Surface Pres- 
sure - Standard lase 

= -?--H 4 shows the contour c~ a .jw 

3i.x waves, the dispLaca.tant body and fte‘“ 
viscous and i.tviacid pressure distributions 
” - —8 standard flow conditions. The ii"- 
-erance i.t t.ta i.tviscid i,td viscous pressures 
Siicws iha jzrsc- ind r.eed *3^ 
interaction. The pressure is ta_culated iron 
in ipprox33ia33cn zo zns ? 3 md 3 i.*“'*^sv* t — 

3ion , iccuxate to aecend order m flrw 
ci.^ation ingia. The inviscid pressure is 
-a«_ulawad isi.tg t.ta local oody slope, wr.ereas 
-^•8 /iscous pressure is oocai.ted cy usi.tg 
slope of t.ta tisplacam.ent body a * 

"'•j *■ ■ "^-e iicferenca i.t the viscous and 

1-. /iscid pressure is due to the differa.nca i.t 
amplitudes of me actual and iisplaca.tent 
occy . .t IS i.ttarasti.tg to coserve t.tat tne 
X1SCOU3 pressure is al.tost teriodic even 
thoug.t t.ta a'/erage displaca.m.ent ftic.tness 
cacraasas . Pigure 5 shows witn the iistri- 
oution of pressure the corresponding listr* - 
= '^-4.on of surface .teat transfer and sxZ.t 
tricticn at tne same base flow bcnditibr.s. 

The pressure peaxs and pea.<s in haati.tg occur 
at about tne sate Location anead of t.te coev 
surface peak. The peak in skin friction is' 
^4icted I.t tne oocosita direction. ;vh • ’ ^ 
the pressure distriouticn is tearlv te-'o-"* - 
tne heating le'/eis and tne skin friction 
pea.ks rise i.t t.ta dewnstraat direction. The 
rata of rise in peak heating is fecreasinc' ' 
vec-y slowly. These results are in con--a2 
-istincticn to our si.tilar studv (3) 3 -« -f- 
.a.minar rlow over a tram of s ine-wav-s , 
wr.ere the peaks i.t heati.tg decreased -->o'd’/ 
i.t the 3traa.trwise direction. rifure 5 ’toiZts 

of t.tk 

3-«.aca parameters is al.tcst u.taf;ac“ad bv 
the presence of additional downstrea.t iistur- 
ba.tce for turbulent boundar/ layers (bbm.oara 


•? 

L 




ilao 5 ind 6). 3ix?iy, Icwr.straam 

wavas iiave ipscraaa iniluar.ee and tne 

prodian saenia Iccaliued. rlaanin? levala ait 
of waves jrow as tl*.a nususer of waves in- 
creases . out t.ta iownstrsajtt si^un friction iS 
unaifectad ov t.ta .vurier of waves. 


\ 

1 — 


3- 


/ 


2i- 

\ 

•Q 

1 — 


01 - 


'*>ji « 153 

iC 32 « iC*"^ 
• aai 


3»^» 138C.TI 
97 «jrra / 


-lot 9 ' 0 t» 

— “ 3i» Waves 



?lg. 5 Sorraca Prccairtias - Standard Oase 



5e.3-lC.32 t iC"'Ti 



sj 0 ^ 


NO ^ove — 

Ore Wove ^ 

^ 'V 

— ^0 A^av<?3 





“ic. 3 ^cvr.straan Sfiacts cn Surraca 
Properties 

To dsr.cnstrata tns affact cn aurraca 
crooartias due to 'dacn -.utaar, vail tanpera- 
fure , and .Peyr.olda nr-tne r , tnrae additional 
cases are snown in Tipures T t.troupn J. The 
i.torsase of Macn .-.’unner (ric’ira 7) crorti 2.5 
to J.5 causes a decrease in tno ratio or 

n .h- . .5 l 3 in t-te standard flow rase, 

tax : . 0 . 

the location vnic.t t.te first wave was o:.aced 
was cheser. in such a way that tne rlat plate 
ooundarv layer displacatent thickness was 
aoout the 5a.te as in t.ts co.T.parison experi- 
■te.ttal study of Reference (^) . 


V' 


The Icwerinp of wail tetoerature wO 
» 0.40 (rip’ure 3) shows a sini.'.ar trend 


"ta.'C "f . 0 . 


IS for the i.tcrease in 


n'urJcer. But the aosoiuta rata of surface 
heating is .tuch hig.ter than i.t the previous 
case. I.tterestingiy , t.te h.'hj^^ curve is 
sr.oother hers than in ether rasas . 

lastly, an increase in Reynolds r.'umer, 
shown in "igure 3, is seen to cause an i.t- 
creasa in t.te ratio of oeux heaci.tc. 
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?ig. T hacn Nuxeer Tffact cn Surface 
Properties 
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?ig. 3 wall Tar.oeratura Effect cr. Surface 
Properties 
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rig. 3 Pe'/nolds Numer 2: feet cn Surface 
Procerties 
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An uiterescincf is?«cr oi tiia present 
results 13 tne Itcacion tf tne peaxs in 
prsssurs, neat trmsiar and shear. The pre- 
sent predictions show the peax3 in pressure 
and neating to cccu.r at asout tae sane loca- 
tion. This 13 in agreement with e.xperimentai 
oo3ar*/ations (^) . The location of the peax 
pressure i.n t.he present results is s.hiitad 
to the right of the location of the i.nviscid 
?e£ix pressure location ( at / • 3) oy a phase 
angle of aaout 50®. Theoretical studies ay 
Inger and Williijivs (^) and ay lexoudia at 
al. (^) predict such a sniit. Data from 
these studies given up to » 2.3 snow a 
shift to the .eft wmoh drops off cuio.*cly 
towards tero at » 2. It is therefore 
?cs*io_e to expect a pnase angle i.n the 
opposite sense for > 2, as is the rase in 
present results . The tiaximua vail shear 
location obtained from present calculations 
IS sniftad to the right of the peaxs in 
pressure and surface neati.ng oy ahcut 50®. 
•Kccordi.-.g to theoretical predictions (13} 
cuA-i-itati-.'eiy suo.n a shift is e-xpectad. 
Txperir.a.ntal fata availaaie at t.ne sane flow 
conditions (1^) snow a periodic trend i.n 
surface orsssure as veil as _n the surface 
heati-ng iistriouticn . Ths periodic trend in 
surface pressure is ooser’/’ed also in the pre- 
sent predictions virh peax vaiies of 
o p ^ _ » 1.3 at M a 2.5 


C0NC1U3ICNS 

\ nuaerical met-hod capable of nandiing 
•tuitiple intaracti.ng flow regions was 
adapted to tne proo.a.t of t.hic.x turbulent 
boundarv layer over a wa'/v wail. 

The results of oalouiations presented in 
terms of .surface pressure, sxin iriotion , and 
heat transfer distributions disclose features 
distinotly iiffarsnt froa the lasinar case. 
The present results shew a snift in the 
location of t.-.e viscous pressure peaxs rela- 
tive to the peaxs in tne mviscid pressure 
and to tne pe^s in tne wall shear. These 
pnase shifts axe in cuaiitativa agreement 
with theoretical oracictions oasec or. small 


le .ocation of pea.xs in 

fa r 1 rf? • f e 


disturbance t.heor 7 . 

'lacous pressure and .teat transfar coi.ncida, 
and the longitudinal pressure variation is 
periodic. This is i.t acree.ment wit.-, tne 


e.xba rimer. tal data. 

?e nodi Cl tv in aurdaca 


The experiments a.so shew 


.taating distrimuticn , 


• max' - c . o . 
L.2 at A 


and o /o. 

■ max' - : . o . 


3.5. 


.farerence 


:il) the 


ceax 


’■'ni-e the prese.nt res'ults predict a ccr.ti.tu- 
ous increase in .taating indicating a cossio-a 
'w«a.xr.es3 in t.he firmulence medal for a 
surface 'witn raaidl'/ '/erving cur',*et'ure . 

AaCICWTHCGE:<S:iT3 

This wor;< was s'upcortad oy ;iAjA mder 
3rant (ISG 1201 witn hr. lames Ounavi-t of 


While the .taating 


values for oressure are ci'.-e.''. on_v at 
-=max'?f.o. 

distributions in the experiments of P.aferance 
( 11) are .tearl'/ repetitive over cor.sacuti'/e 
waves ('With 


the pangiay .=,esearch 


ntar as 


rac.tr.icaJ 


. s 

.max' ■ f . ? . 

1.2 at k 


‘ 1.3 at *2.5 and 

■ 

3.5) tne present 
.tuous i-tcraasa ever 


■'max' "f . o . 

predictions snow a con 
the length of tne waves. Mote tnougn that 
in. the expenme.ttal study there is also 
indication that saparatio.t cccurs , while a 
Lac.x of separation _s oosarved i.t the anal'/- 
tical res'ults cf figure 5. The ca'usa of t.tis 
disagrae.'T.ar.t is tot certai.t out it could, well 
be due to cur choice cf t'lrbulance mo cal' 
cr in tne fact tnat the cresent 
da tot simulate we .3. ancagn all 
conditions ( t.hree-dimar.siona_ e 
dO’indajT’ Layer development on t 
wall) . Mote t.-.e calculated bou 
dispiacament thic.tn.ess .cf the 
file at station x » '1.30 is 2. 
to the '/a-ue gi''er. in Peference 
aver tne prediotad surface .teat 
this station is too hign when c 
experi.mantal data of .fefare.mce 
dieted value is .tf.g. • 35.5 -wa 
5 2.5 'watts /m- ®t< i.t' experiments 
diction .met.mods also typically 
to about t.he same level tne .teat transfer 
rates for ooundac"/ Layers devel; 
t.he wind funnel walls (^) thus 
fiat some fi.-.al acD'ustm.e.mts m.ay 
t.te furbule.mca .model far these flows. 
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